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Abstract 

The mean kinetic energy of protons in water is determined by Deep Inelastic Neutron Scattering 
experiments, performed above and below the temperature of maximum density and in the super- 
cooled phase. The temperature dependence of this energy shows an anomalous behavior, as it 
occurs for many water properties. In particular two regions of maximum kinetic energy are identi- 
fied: the first one, in the supercooled phase in the range 269 K - 272 K, and a second one above 273 
K. In both these regions the measured proton kinetic energy exceedes the theoretical prediction 
based on a semi-classical model. Noteworthy, the proton mean kinetic energy has a maximum at 
277 K, the temperature of the maximum density of water. In the supercooled metastable phase 
the measured mean kinetic energy and the proton momentum distribution clearly indicate proton 
derealization between two H-bonded oxygens. 

PACS numbers: 61.25.Em,78.30.Cp,78.70.Nx 



[*] Author to whom correspondence should be addressed. E-mail: bruni@fis.uniroma3.it 



1 



Among liquids, water is undoubtely the most studied, as it is the medium of life and 
chemistry on our planet, and its puzzling properties are a challenge for science Q, 2]. Yet 
water quantum behavior has received less attention than its classical thermodynamic prop- 
erties, in spite of the modest cooling needed to evidence quantum effects in water. These 
have been recently suggested to determine the water density maximum 3j and shown to 
influence the molecular geometry Q] and in particular the OH bond-length. Quantum ef- 
fects have also been evidenced in the short time dynamics of protons along the OH bond in 

by Deep Inelastic Neutron Scattering (DINS) [s| 



supercooled 5] and confined water 
experiments and computer simulations [10 



111 ]. Here we report DINS measurements on bulk 



water in the stable phase around the density maximum and in the supercooled metastable 
regime, showing clear evidence for quantum effects on the momentum distribution n(p) and 
mean kinetic energy (Ek) of water protons. 

At present DINS is the only technique which determines the momentum distribution and 
mean kinetic energy of a single particle in condensed matter. It is based on neutron scattering 
measurements at high energy (1 eV < hu < 200 eV) and momentum (120A -1 < hq < 
300A _1 ) transfers, thus probing both the short-time dynamics and local environment of the 
particle. The high energy and wave vector transfers achieved allow the scattering process to 
be described within the framework of the Impulse Approximations (IA) [12.]. This assumes 
that neutron scattering is incoherent and occurring within time scales much shorter than 
the typical relaxation times of the collective excitations of the system. In these conditions 
the struck particle recoils freely from the collision, with inter-particle interaction in the final 
state being negligible (i.e. the wave-function of the particle in its final state assumed to be a 
plane wave). In a molecular system, as for instance water, the contribution to the scattering 
cross section of protons can be easily distinguished from that of oxygens, due to the large 
mass difference. The scattering cross section is then expressed in terms of the single proton 
momentum distribution n(p), whose variance is directely related to its mean kinetic energy 
(Ek). These quantities provide a richness of information about the potential surface that 



the proton experiences |13j , inc 



uding 



14 



the effects of hydrogen bonding, thus complementing 



15] and allowing a direct comparison with quantum 



microscopic structural studies, 
Monte Carlo simulations [h], fill ]. 

DINS experiments have been carried out on the VESUVIO spectrometer [9], at the ISIS 
spallation neutron source (UK). The samples were contained in the same disk shaped alu- 
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minum can (5 cm diameter, 1 mm thickness) with inner Teflon coating used in previous work 



n 



on supercooled water [5|. Experimental data have been recorded at 64 fixed-angle detectors 
in the standard Resonance Detector configuration 16j employing the Foil Cycling technique 



171 ] . Data have been corrected for multiple scattering, aluminium and oxygen signal, by 



using the same routines as in previous work [5|. Within the IA framework the dynamical 
structure factor, Sia(q,^), is related to the n(p) through the relation: 



S IA (q,uj) 



n(p) 8 I uj 
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dp 
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where tko is the energy transfer, is the recoil energy of the struck atom of mass M, and q 
is the vavevector transfer. The standard deviation of the n(p) function, a, is related to (Ek), 
through the relation a 2 = |jjf (E^). The dynamical structure factor is then expressed in terms 
of a Neutron Compton Profile (NCP) fl: F(y) = ^S IA (q, u), where y = g (u - gj) is the 
West scaling variable js]. The NCP lineshape is convolved with the instrumental resolution 
function, and represented as a series expansion in Hermite polynomials: 
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The coefficients c n and a, appearing in the series expansion, have been determined by a least 
squares fitting procedure (see caption of Fig. 2), and small corrections due to deviations 
from the IA were taken into account [9fl. We did not attempt to correct for the so-called 

nnn 

intensity deficit seen for hydrogen relative to heavier nuclei [18|, ]19|, ]20( ■ Altough the issue of 



y 



(0.2) 



the intensity deficit in DINS experiment is still debated, it should be noted that in systems 
where the n(p) and the NCP are purely gaussian, such as ZrH 2 , no distortion of the measured 



NCP is indeed observed 



2l|. Finally the n(p) is expressed by the expansion 
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where Ln are generalized Laguerre polynomials. Full details about data analysis may be 
found in ref. ]9|. 

Figure 1 shows the temperature dependence of proton's (Ek) measured in both stable and 
metastable phases of bulk water in the temperature range from 300 K to 269 K, along with 
the semi-classical prediction, E sc . The temperature dependence of E sc in the stable water 
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phase, shown in Fig. 1 as a dotted line, has been derived in Ref. |22|, taking into account 



translational, rotational, and vibrational contributions, including the optical spectroscopic 
data available in the literature. While data for ice 23| and water [16j above 293 K, including 



supercritical states 



22 



241 ]. are satisfactorily described by the semi-classical prediction, data 



around the temperature of maximum density and in the supercooled phase show an excess 
of proton mean kinetic energy. In particular (Ek) shows two maxima: one at 277 K and 
the other in the supercooled phase in the range 269 K - 271 K, with an excess of (Ek) 
with respect to E sc of about 3 kJ/mol (30 meV) and 12 kJ/mol (120 meV), respectively. 
The peculiar temperature dependence of (Ek), suggests that two distinct mechanisms may 
be considered below and above the melting temperature. As a matter of fact the excess 
of (Ek) in the stable water phase is moderate and its temperature dependence follows that 



of density 



251 ]. showing a maximum at the same temperature (see insert in Fig. 1). A 



correlation between density and (Ek), as measured by DINS, is not a novelty in principle, 



as it has already been observed in helium [26|, [27[ and explained in that case by using a 



harmonic model for the fluid. Nevertheless the peculiar density behavior of water and the 
temperature evolution of the H-bond network does not allow a similar model to work in this 
case. We suggest instead that the anomaly of (Ek) vs T above 273 K may be explained 
as a further evidence for water structural anomalies, which manifest through the existence 



of a maximum of density and transport properties in the stable water phase 



28|. Within 



this hypothesis the maximum of (Ek) at 277.15 K, shown in Fig. 1, may be an indirect 
manifestation of the competition between zero point energy, Eq, and thermal fluctuations, 
which has been proposed as the quantum origin of the density maximum [3J]. We notice 
however that quantum effects are not necessarily required to explain the existence of a 



maximum of density in water 



29J]. 



On the other side the huge increase of (Ek) in the metastable states of water, already 
observed in a recent publication 5] and confirmed by present new measurements at 272.15 
and 272.95 K, seems to be directely related to the likely derealization of protons along the 
H-bond. Fig. 2 shows the radial momentum distribution, 4irp 2 n(p), at 277.15 and 271.15 
K. This function at the lowest temperature shows a narrowing at \ow-p and a shoulder at 
high-p (~17A _1 ), compared to that measured at 277.15 K, compatible with the transition 
from a single to a double well potential. The presence of a shoulder at high-p in the radial 
momentum distribution has indeed been ascribed to coherent derealization of protons over 
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two sites of a double well potential 13] felt by the proton along the H-bond direction between 
two water molecules. In other words, a statistical ensamble of water molecules in which 
protons are localized in the vicinity of the covalently bonded oxygens is predicted to show a 
single maximum of Anp 2 n(p) at p ~ 6 A -1 , corresponding to the intramolecular O-H distance 
d = 2j£ ~l A, as shown in Fig. 2 and confirmed by quantum mechanical simulations id. 111]. 
Conversely, the appearence of a shoulder at p ~ 17 A -1 in the metastable supercooled 
phase indicates a derealization of the proton over a distance Ad = ^ ~ 0.6 A from the 
equilibrium position. We notice that (E k ) in the supercooled phase is comparable with the 
H-bond energy (~ 20 kJ/mol) and Ad is ompatible with the width of the fluctuations of 
the H-bond length (width of the first intermolecular peak of the oxygen- hyd rog en radial 
distribution function) and the oxygen-oxygen distance in supercooled water 14]. Based 



on this observation, we proposed 



5] the quantum excess of proton mean kinetic energy 



be correlated to the average distance between two first neighboring H-bonded oxy 



ens m 



water. Such assignement was corroborated by similar evidences in confined water 



15|- 



Further evidences for a correlation between proton derealization and distance between its 



neighboring oxygens are given in Ref. 



31 



showing the changes of the potential 



3d and 

landscape, energy barrier and wave-function of a proton as a function of this 0-0 distance. 

The present results show a clear and unexpected anomalous temperature dependence of 
the proton mean kinetic energy. Below 273 K, in the supercooled phase, this anomalous 
behavior is associated with a coherent derealization of the proton between first neighbor 
oxygens. Above 273 K the temperature dependence of the proton mean kinetic energy 
resembles that of water density. 
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Figure Captions 



Figure 1: Water proton mean kinetic energy (Ek) as a function of temperature. Open 
symbols have been used for water in the stable phases, solid symbols for the metastable 



supercooled phase, namely: ice 



at 269 K (square); liquid stable water (circles). 



16 



those below 272.15 K are from Ref. 



Water data above 293 K are taken from Ref. 

[5I. Present data are reported for temperatures in the range 272.15 - 285.15 K. The 
dotted line corresponds to the semi-classical prediction [3] . The error bars are derived 
from the least squares fitting procedure used to derive (E K ) from the differential cross 
sections (see references (| M). The solid and dashed lines are guides for the eye. In 
the insert (E K ) data in the stable liquid phase (circles and right axis) are reported 
in comparison with the density of water (line plus solid triangles and left axis) as a 
function of temperature. 

Figure 2: Spherically averaged momentum distribution of water protons at the point 
of maximum density (present experiment, T = 277.15 K), (dotted line) and in the 
metastable supercooled phase at T = 271 K (solid line, from Ref. |5|). Experimental 
uncertainties are less than ±1 %. Values for the fitting parameters (see Eq. 2-3) are 
as follows: C\ is set to zero by definition, the non-Gaussian coefficient, C2, and a are 
(0.443 ± 0.008) and (6.05 ± COS)!" 1 at T = 271 K; (0.148 ± 0.010) and (5.13 ± 
0.02)A _1 at T = 277.15 K; the coefficients c n > 3 are found to be negligible. 
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